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Random walks (RW) of particles adsorbed in the internal walls of porous deposits produced by 
ballistic-type growth models are studied. The particles start at the external surface of the deposits 
and enter their pores, in order to simulate an external flux of a species towards a porous solid. 

For short times, the walker concentration decays as a stretched exponential of the depth z, but a 
crossover to long time normal diffusion is observed in most samples. The anomalous concentration 
profile remains at long times in very porous solids if the walker steps are restricted to nearest 
neighbors and is accompanied with subdiffusion features. These findings are correlated with a decay 
of the explored area with 2 . The study of RW of tracer particles left at the internal part of the solid 
rules out an interpretation by diffusion equations with position-dependent coefficients. A model of 
RW in a tube of decreasing cross section explains those results by showing long crossovers from an 
effective subdiffusion regime to an asymptotic normal diffusion. The crossover position and density 
are analytically calculated for a tube with area decreasing exponentially with 2 and show good 
agreement with numerical data. The anomalous decay of the concentration profile is interpreted as 
a templating effect of the tube shape on the total number of diffusing particles at each depth, while 
the volumetric concentration in the actually explored porous region may not have significant decay. 

These results may explain the anomalous diffusion of metal atoms in porous deposits observed in 
recent works. They also confirm the difficulty in interpreting experimental or computational data 
on anomalous transport reported in recent works, particularly if only the concentration profiles are 
measured. 

PACS numbers: 05.40.-a, 66.30.Pa, 68.55.-a, 81.05.Rm 

I. INTRODUCTION 

Diffusion in porous media long has been a topic of in¬ 
terest [ 3,0 due to a variety of technological applications, 
which are usually related to the large surface-to-volume 
ratio of those materials. The concentration profile of a 
diffusing species entering a disordered medium is typi¬ 
cally of the form 

p( 2 ,t)-exp[- 2 '"/R(t)], (1) 

where z is the direction perpendicular to the external 
surface, with z = 0 at that surface. The function R (t) 
measures the spreading of the concentration profile in 
time t and asymptotically behaves as 

( 2 ) 

From Eqs. © and the mean-square displacement of 
the diffusing material scales as 

(z2) _ (3) 

Normal (Fickean) diffusion is characterized by a spatially 
Gaussian distribution of the diffusing species and linearly 
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increasing mean-square displacement: /i = 2 and a = 1. 
Otherwise, the diffusion is called anomalous. The scaling 
of (z^) is frequently used to separate cases of subdiffusion 
{a/^ < 1/2) and superdiffusion (a/@ > 1/2). The most 
frequently observed anomaly is subdiffusion because the 
irregularities of the porous media (e. g. barriers and dead 
ends) restrict the movement of the diffusing species. 

The relevance of anomalous diffusion in real systems 
has been shown in several recent works. For instance, 
fluorescence spectroscopy was used to distinguish cases 
of normal and subdiffusion of tracer particles in colloidal 
crystals 0, and combined with confocal microscopy (to 
probe concentration profiles), it was used to show that 
reaction-diffusion models were necessary to explain dye 
transport in metal-organic frameworks 0. Nuclear mag¬ 
netic resonance was recently used to measure exponents 
a and @ of water diffusion in several disordered col¬ 
loidal systems, showing deviations from Fickean diffusion 
0 0 . Particularly interesting cases of stretched exponen¬ 
tial concentration profiles and superdiffusion were shown 
when Pt atoms entered the pores of porous carbon and 
anodic porous alumina during plasma sputtering deposi¬ 
tion on those samples 00. Moreover, due to the impor¬ 
tance of subdiffusion phenomena in cells and model sys¬ 
tems, the need for experimental standards was recently 
highlighted M- 

There is also an intense theoretical work on anomalous 
diffusion. Some approaches involve the computational 
study of fluid flow and solute transport in disordered 
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media [HHii , analytical methods based on advection- 
diffusion equations [id . fl^ and fractional diffusion equa¬ 
tions [3 . The study of random walks (RW) in 

models of disordered media Qi is another simple and 
widely used approach, whose relevance to describe real 
system features is illustrated in recent works 

Motivated by the observation of anomalous diffusion of 
Pt atoms penetrating samples of disordered porous car¬ 
bon 0, in the first part of this work we study RW of 
particles adsorbed in the internal surface of porous de¬ 
posits after being released at their outer surface. The de¬ 
posits are produced by ballistic deposition (BD) [H, 
and by an extension of that model [281 - 1^ . For short 
times, the concentration profiles are stretched exponen¬ 
tials [slower than Gaussian, jj, < 2 in Eq. ([T])]. This fea¬ 
ture is still observed at long times in one type of deposit, 
and the corresponding time scaling of the concentration 
profile is consistent with subdiffusion. The area explored 
by the walkers decreased approximately as the inverse of 
the depth in this case. Tracer particles inside the samples 
do not show dependence of diffusion coefficients on the 
depth, which rules out theoretical approaches with this 
assumption. Although the values of exponents fi and a 
differ from those of Ref. Q , the scaling of the concentra¬ 
tion with depth and time show similar deviations from 
normal diffusion. 

The second part of this work is devoted to the study of 
random walks in tubes of decreasing cross section, which 
show concentration profiles with stretched exponentials 
and anomalous scaling of the mean-square displacement. 
A tube with area decreasing with the inverse of the depth 
(similarly to the porous deposits with anomalous diffu¬ 
sion) shows significant deviations from normal scaling 
and exponent values much closer to those of Pt atoms 
entering porous carbon samples Q- A detailed discus¬ 
sion of the crossover from anomalous to normal diffusion 
is presented for the case of tubes with exponentially de¬ 
creasing cross section, in which analytical solution is pos¬ 
sible. These results confirm the difficulty in interpreting 
experimental or computational data on anomalous trans¬ 
port, particularly the concentration profiles, consistently 
with recent works [3, [l^, [H, [s^ ■ 

This paper is organized as follows. In Sec. II, we 
present the models of ballistic deposits and the random 
walk models. In Sec. Ill, we study the scaling of concen¬ 
tration profiles inside the ballistic deposits and related 
quantities. In Sec. IV, we study diffusion in tubes of 
decreasing cross section and discuss the crossover in con¬ 
centration scaling. In Sec. V, our conclusions are pre¬ 
sented. 


II. POROUS DEPOSITS AND RANDOM WALK 
SIMULATION 


Here we describe the two types of disordered porous 
media used in this work and the two types of RW simu¬ 
lated inside those samples. 



FIG. 1: (Color online) Illustration of the aggregation rules 
and cross sections of the tops of deposits generated by (a),(b) 
BD and (c),(d) BDNNN. Dark gray cubes are previously ag¬ 
gregated solid particles and light gray cubes are incident par¬ 
ticles, whose aggregation positions (also shown in light gray) 
are indicated by arrows. 


In the BD model, particles are released from a ran¬ 
domly chosen position above a d-dimensional substrate, 
follow trajectories perpendicular to an initially flat sub¬ 
strate [here a (x, y) plane] and stick upon the first con¬ 
tact with a nearest neighbor (NN) occupied site, which 
may be the substrate or a previously deposited particle 
[^[^. The aggregation rules are illustrated in Fig. la. 
The resulting aggregate is porous and has a rough sur¬ 
face, as illustrated in the cross-sectional view of a three- 
dimensional deposit in Fig. lb. Long time simulations 
indicate that the porosity is approximately 0.667. 

A simple extension of that model, called BDNNN, al¬ 
lows aggregation of the incident particle by contact with 
a NN or a next nearest neighbor (NNN) occupied site 
[2§ - l3ll| . as illustrated in Fig. Ic. The resulting deposit 
has a larger porosity, approximately 0.834, as illustrated 
by the cross-sectional view of Fig. Id with an identical 
number of deposited layers as in Fig. Ic. 

Deposits of average heights near 900 lattice units and 
lateral sizes 1024 were grown. The structure of each sam¬ 
ple remains fixed during the simulation of diffusion, i.e. 
they are nondeformable porous solids. A single deposit 
is large enough to represent all the microscopic environ¬ 
ments that are relevant for the RW statistics. For this 
reason, only 3 samples were produced by each growth 
model (BD and BDNNN) and used in RW simulations, 
providing approximately the same average quantities; for 
instance, the dispersion in the porosity is smaller than 
0.3%. 

The maximal height of a solid particle in each column 
(x, y) of a deposit is defined as the column height h(x, y). 
The set {h{x,y)} defines the external surface of the de¬ 
posit. The average h of that set is taken as the position 
z = 0 for RW simulations. The z axis is oriented to the 
interior of the deposit, perpendicularly to the substrate 
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FIG. 2: (Color online) Illustration of possible steps of RW 
of the species (dark blue) moving inside the porous medium 
(solid part in dark gray) according to rules (a) (I) and (b) 
(II). The target position of the allowed steps is shown in light 
blue. 


where it was grown. This is illustrated in Figs, lb and 
Id. 

In each porous sample, 2^® (~ 10®) non-interacting 
walkers are left at t = 0 at randomly chosen points with 
z = 0, with the condition that at least one NN of the 
starting point is a solid site. In one time unit, one step 
trial is performed for each walker. Adsorption to the 
porous solid is always required, i. e. only steps to points 
that also have a NN solid site are allowed. 

Independent simulations were performed with two pos¬ 
sible conditions for choosing the steps of the walkers: (I) 
a step to a NN site of the underlying cubic lattice (6 pos¬ 
sibilities) is randomly chosen and (II) a step to a NN or a 
NNN site of the underlying cubic lattice (18 possibilities) 
is randomly chosen. If the target site has at least one NN 
of the porous solid (adsorption condition), then the step 
is executed. Otherwise, the step trial is rejected and the 
walker remains at the same position. 

Figs. 2a and 2b illustrate steps with conditions (I) 
and (II), respectively. Diffusion under condition (I) is 
severely restricted because it does not allow corner round¬ 
ing. However, with condition (II), this process is possible 
and steps from one branch of the porous solid to another 
one are facilitated. 

The maximum simulation time is 10®, corresponding 
to the maximum number of possible steps of each walker 
(due to the no-desorption condition, the number of exe¬ 
cuted steps may be much smaller than that). This maxi¬ 
mal time is suitable to avoid that any walker reaches the 
bottom of the substrate. 

As a final remark, we recall that there are other ex¬ 
tensions of BD to represent a variety of porous materials 
(33l - l^ . However, the aim of this work is to understand 
basic features of anomalous diffusion of a species enter¬ 
ing a porous material, instead of representing a particular 
application. For this reason, we restrict our work to the 
study of samples produced by BD and BDNNN models. 


III. DIFFUSION IN THE POROUS DEPOSITS 
A. Scaling of concentration profiles 

Figs. 3a and 4a show log [—log (p/po)] as a function 
of log 2 in three different times for RW in porous solids 
produced by BD (lower porosity - Fig. lb), respectively 
with conditions (I) and (H) for the steps. Here, po is the 
maximum value of the density profile. Figs. 5a and 6a 
show the same quantities for RW in porous solids pro¬ 
duced by BDNNN (high porosity - Fig. Id), respectively 
with conditions (I) and (H) for the steps. 

Linear fits of the data for large z are also shown in 
Figs. 3a, 4a, 5a, and 6a. Their slopes are estimates of 
the exponent p in Eq. [T] At short times (t = 10®), 
the slopes are significantly below 2 in all cases. This is 
usually a transient behavior, since the fits of the data at 
longer times (t = 10"^ and t = 10®) in Figs. 3a, 4a, and 
6a have slopes close to the normal diffusion value p = 2 
(deviations are smaller than 10%). However, in BDNNN 
deposits, the stretched concentration decay (slower than 
Gaussian) is observed up to very long times, with fits in 
Fig. 5a ranging from p Ri 1.33 to p ~ 1.40, in depths 
ranging from 2 ~ 1 to z Ri 250 (in units of the lattice 
constant). 

The plots of the concentration profiles are rescaled by 
the time t in Figs. 3b, 4b, 5b, and 6b. Reasonable data 
collapse of the large z data is obtained in Figs. 3b, 4b, 
and 6b, using the scaling variable which gives 

a/p = 1/2 [Eqs. (P) and (P)]. This is consistent with 
asymptotic normal diffusion. 

On the other hand, walkers of type (I) in BDNNN de¬ 
posits show much slower time dependence of the concen¬ 
tration profile, which could be antecipated by comparison 
of Fig. 5a with Figs. 3a, 4a, or 6a. A reasonable data 
collapse for large z is obtained in Fig. 5b with the scal¬ 
ing variable z/t® ®®, which gives a/p ~ 0.15 [Eqs. ([T]) 
and ©]■ This estimate gives (z^) ^ t® ®, which indicates 
subdiffusion. 

A comparison with experimental results is interesting 
at this point. 

For Pt atoms entering porous carbon samples. Ref. 
0 gives p Ri 0.55, which characterizes a stretched decay 
(much slower than Gaussian) of the concentration profile, 
and gives a very slowly increasing R{t), with a r; 0.2. 
Gombination of these results gives (z^) ~ i. e. 

0.36 [Eqs. (P and (P)]. In the usual classification 
Ul 0], a system with a/p < 1/2 is a case of subdiffusion 
(however, the remarkable stretching of the concentration 
profile is a nontrivial feature that led the authors of Ref. 
[7| to propose it was superdiffusive scaling). 

For water diffusion in disordered colloidal systems, Ref. 
0 shows that several combinations of a < 1 and p < 2 
are possible, which may give sub or superdiffusion. How¬ 
ever, in that case, all estimates are close to a = 1 and 

p = 2. 

A recent work on Pt atoms entering pores of anodic 
alumina 0 reported p Ri 1/3 and a Ri 1.25, which gives 
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FIG. 3: (Color online) Concentration profiles without (a) and 
with (b) time scaling, number of sites explored by RW (c), and 
local diffusion coefficients (d) of tracer particles as a function 
of the depth z in porous solids produced by BD with rule 
(I) for the random steps of the RW. In (a), (b) and (c), red 
squares, green crosses and blue triangles relate respectively to 
t = 1000,10000,100000. Full, dashed, and dash-dot-dot lines 
in (a) are linear fits for the largest values of z. 


FIG. 4: (Color online) Concentration profiles without (a) and 
with (b) time scaling, number of sites explored by RW (c), and 
local diffusion coefficients (d) of tracer particles as a function 
of the depth z in porous solids produced by BD with rule (II) 
for the random steps of the RW. Points and lines are defined 
in the caption to Fig. [3] 
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FIG. 5: (Color online) Concentration profiles without (a) and 
with (b) time scaling, number of sites explored by RW (c), and 
local diffusion coefficients (d) of tracer particles as a function 
of the depth 2 in porous solids produced by BDNNN with 
rule (I) for the random steps of the RW. Points and lines 
have identical meaning as in Figure |3l 


FIG. 6: (Color online) Concentration profiles without (a) and 
with (b) time scaling, number of sites explored by RW (c), and 
local diffusion coefficients (d) of tracer particles as a function 
of the depth 2 in porous solids produced by BDNNN with 
rule (II) for the random steps of the RW. Points and lines are 
defined in the caption to Fig. [3] 
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a remarkably rapid superdiffusion [2a//i Ri 7.5; Eq. Q]. 
However, the organized geometry of those samples make 
them differ markedly from our deposits. 

Finally, we recall that subdiffusion was also observed 
in porous deposits produced by other ballistic-like models 
in Ref. [l^. Those deposits seem to have lower porosity 
than ours and have different pore shapes and connec¬ 
tivity. Only the mean-square displacement of the walk¬ 
ers was studied in that work, thus the present results 
strongly suggest the study of concentration profiles in 
those samples. 


B. Explored area and diffusion coefficients 

The anomalous diffusion for walkers of type (I) in 
BDNNN deposits is related to several constraints for their 
steps, similarly to other systems 0. First, the volume 
accessible for the walkers is restricted because the poros¬ 
ity is very large (Fig. Id) and they are always adsorbed 
to the internal solid walls. Second, since the steps are 
restricted to NN (Fig. 2a), some movements become im¬ 
possible, such as the corner rounding shown for the upper 
walker in Fig. 2b. This is certainly an important contri¬ 
bution to anomalous scaling. 

Here we measure some quantities that help to explain 
the results presented in Sec. IHI Al and motivate the mod¬ 
els presented in Sec. El 

The pore region accessible for the walkers is character¬ 
ized by the total area A{z) explored by the walkers at 
each depth z, for several times. This area is the total 
number of pore sites at depth ^ that have been occupied 
by a walker at least once up to time t. The explored areas 
A{z) for three different times are shown in Figs. 3c, 4c, 
5c, and 6c, for the respective solids and step conditions. 

The three cases with transient anomalous scaling (Figs. 
3, 4, and 6) show similar behavior of A{z) for short and 
long times. For short times (t % 1000), there is a rapid 
decay of A{z) at depths varying from r; 40 to ~ 100, de¬ 
pending on the type of solid and the allowed steps. The 
corresponding fits of the concentration profiles extend to 
larger depths and show anomalous scaling (p, < 2). For 
long times, the flat region of A{z) extends to much larger 
depths, which means that the walkers explore an approx¬ 
imately constant cross section as they penetrate in the 
porous deposits. The corresponding concentration pro¬ 
files show normal scaling fj. = 2. 

For walkers of type (I) in BDNNN deposits. Fig. 5c 
shows a slow decay of A{z) until very long times. A fit of 
the data in Fig. 5c for t = 10^ suggests an approximate 
decay as A{z) \/z. Thus, there is an actual reduction 
of the area that can be reached by the walkers as they 
penetrate in the porous solid. The porosity is approxi¬ 
mately depth-independent, which means that a large part 
of the porous space is not accessible to particles moving 
with the constraint of adsorption to the internal walls. 

Figs. 7a-d show cross sections of a BDNNN solid at 
four depths z and highlights the area occupied by walk- 
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FIG. 7: (Color online) Cross sections of a porous solid pro¬ 
duced by BDNNN at depths (a) « = 16, (b) 32, (c) 96, and 
(d) 128, with solid part in black, porous sites visited by RW 
[rule (I)] at t = 10® in red, and non-visited porous part in 
white. 


ers of type (I) at t = 10®. It confirms the reduction of 
the explored area with the depth. As will be shown in 
Sec. E the stretched decays of concentration profiles is 
intimately related to this feature. 

Many works on anomalous diffusion connect this fea¬ 
ture to position-dependent diffusion coefficients (see e. 
g. Refs. [H,!!^). In the present systems, this possibil¬ 
ity is ruled out by measuring local diffusion coefficients 
that characterize the diffusion process in a narrow range 
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around each depth z (in opposition to asymptotic coef¬ 
ficients that would characterize the diffusion in the full 
pore network). 

The local diffusion coefficient at height z is measured 
by tracers left at all accessible points in the correspond¬ 
ing depth of the solid. The number of tracers left at each 
point (x, y, z) is proportional to the number of walkers 
which have occupied that position in the course of the 
original simulations up to t = 10^ (i. e. in the simu¬ 
lations with walkers left at the external surface). Each 
tracer tries to execute a small number of steps, to = 100, 
following the same rules [(I) or (II)] of the RW. For the set 
of tracers starting at a given depth z', the mean-square 
vertical displacement ((Az)^) (z',to) is measured after to 
steps (averaging over x, y, and different tracers). The re¬ 
sulting local diffusion coefficient at depth z' is estimated 
as 


^ ^ ,, {iAzf)iz',to) 

D^{z) = --- 

CO 


(4) 


(z') and Dy (z') are equivalently measured. 

Figs. 3d, 4d, 5d, and 6d show (z) as a function of 
the depth z, for the respective solids and walker types of 
Figs. 3-6. Considering the small values of the diffusion 
coefficients, we observe that only a small region around 
the original position is explored, thus we are actually 
measuring local diffusion coefficients. In all cases, the 
coefficients are approximately depth-independent, which 
means that the porous solids are homogeneous at short 
lengthscales. Thus, the anomalous diffusion observed in 
these systems cannot be explained by changes in the local 
diffusion coefficient. This contrasts to some widely used 
models of anomalous diffusion and the proposed 

generalized diffusion equation of Ref. Q ■ 

Figs. 3d, 4d, 5d, and 6d also show the diffusion coeffi¬ 
cients Dx{z) and Dy(z), with no significant dependence 
on z. The estimates of Dx{z) and Dy{z) are nearly the 
same, which confirms the homogeneity in the horizontal 
directions and the accuracy of the simulations. 

We conclude that local features of the porous deposits 
are not responsible for the anomalous diffusion observed 
here. Instead, the anomaly is related to the different 
accessibility of inner sites at a given height, since they 
have different connections (if any) to the outer surface. 


IV. DIFFUSION IN TUBES OF VARIABLE 
CROSS SECTION 


The work with RW entering porous deposits showed 
anomalous diffusion when the available area for the walk¬ 
ers decreased with the depth, but with constant diffusion 
coefficients. For this reason, this second part of the paper 
considers a simple model with those features. 


A(0) 


walkers 
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FIG. 8: (Color online) (a) Tube with z-dependent cross sec¬ 
tion and walkers moving inside it; (b) Random steps of the 
walkers in the tube, with larger number of steps towards the 
larger neighboring areas. 


A. Basic equations and crossover 


We consider an ensemble of RW starting at randomly 
chosen positions with z = 0 at t = 0, confined to move in 
the three-dimensional region z > 0, + y'^ > [A (z) /tt], 

which defines a tube of z-dependent cross section A{z). 
Here we are interested in cases where A{z) is monotoni- 
cally decreasing, similarly to the porous media produced 
by BDNNN. An illustration is provided in Fig. 8a. 

Let N{z,t) be the number of walkers at position z at 
time t. This quantity is proportional to the walker con¬ 
centration p{z, t) if the latter is measured over a constant 
cross section, i. e. the cross section of a solid containing 
the empty tube. This is similar to the measurement of 
concentration in the porous deposits of Sec. IIII Al and in 
experimental works, since the total number of walkers at 
each depth z is divided by the (constant) area of a slice 
that includes solid parts and pores. 

In the limit of continuous z and t, the change in N(z, t) 
in a small time interval is 


N {z A A At) — N (z, t) = 

— kN (z, t) + kN (z — Az, t) 


A{z — 2Az) + A (z) 


kN (z -h Az, t) 


A(z) 

A (z -I- 2Az) -h A (z) ’ 


(5) 


where is a constant (increasing with At and decreas¬ 
ing with Az). Eq. ([S]) accounts for the fact that all 
neighboring sites are equally probable for a random step, 
thus the number of walkers that move from z ± Az to 
a neighboring position (z ± 2Az or z) is proportional to 
the area available at that position. This is illustrated in 
Fig. 8b, with a larger number of walkers moving to the 
larger neighboring area. 

Considering that A (z ± 2Az) -I- A (z) ~ 2A (z ± Az) 
and dividing Eq. ([S]) by A (z), we obtain an equation for 
N (z,t) /A(z). Defining D = kAz^/At and taking the 
continuum limit (Az —>■ 0, At —>■ 0), we obtain 


d{N/A) Dd^N/A) 

dt 2 9z2 ■ ^ ^ 

This is the usual diffusion equation for the density 
N{z,t)/A{z). 











The solution of Eq. m is 




N {z,t) = NoA{z)exY)[-z'^/ {2Dt)], (7) 

where Nq is a normalization constant dependent on t 
(typically as a power law, which is slowly varying com¬ 
pared to the exponential factor); for instance, for con¬ 
stant A, the normal value NnA = , is recovered. 

This solution for N (z,t) is a normal diffusion behavior 
with the concentration at position 2 templated by the 
area of the tube. 

Rewriting the walker concentration of Eq. o as 

Nexp [—z'^/{ 2 Dt) + In A (z)], (8) 

a crossover is expected if | In A (z)| does not scale as i. 
e. non-Gaussian A (z). Here we restrict the discussion to 
the cases in which A (z) decreases slower than a Gaussian. 
In these cases, the first term in the exponential of Eq. dH) 
is dominant for large z, but the second one is larger for 
small z, particularly for long times. A crossover position 
Zc is found by matching the two terms in the exponential 
of Eq. (I5|) , which gives 

ZcV(2^0 «lnA(z,). (9) 

Normal diffusion concentration is observed for z ^ z^ 
and anomalous diffusion is observed for z <C Zc, with 
a stretched exponential decay (slower than Gaussian) of 
the concentration profile. 

As time increases, z^ increases, so the anomalous diffu¬ 
sion regime extends to larger regions and to lower densi¬ 
ties. In this regime, matching the dominant terms of Eqs. 
© and the anomaly exponent is given by power¬ 
counting in the relation In A (z) ^ z^. 

Normal diffusion is able to distribute a large particle 
density in domains of size 

Zdiff ~ VDi, (10) 

thus this is a region far from the Gaussian tail of the 
concentration profile. If diffusion takes place in a phys¬ 
ically confined region, then the particle density (taken 
over porous and solid regions) is templated by the shape 
of the confined pore space up to z ~ Zdiff. In other 
words, the concentration profile scales similarly to A{z) 
in this highly populated region. Gonsequently, it is the 
templating effect of the tube shape that leads to a concen¬ 
tration decay similar to A{z) for small z and the crossover 
phenomenon described above. 

B. Tubes with A{z) ^ Aojz 

As a first numerical application, we consider a tube 
with A{z) = Aq/z, where Aq = 10^. This is a decay 
similar to the explored area of RW of type (I) in BDNNN 
deposits (Sec. IIIII) . Simulations of 10® different walks 



ln(z/t“'‘") 


FIG. 9: (Color online) Normalized walker concentration as a 
function of the scaled depth in a tube with A{z) = Aq/z 
at t = 1000 (red squares), t = 5000 (green crosses), and 
t = 10000 (blue triangles) and the solid line has slope 0.7. 


were performed. In each walk, the particle executes one 
random step to a NN site per time unit, thus H = 1/3 in 

An 

Eq. (|6]). Eq. d?]) gives N ~ — exp [— {^Dt)]. 

In Fig. 9, the usual scaling plot of the concentration 
profile is shown. The scaling variable provides the 

best fit for large z, where the slope is near 0.7, which is 
much smaller than the normal diffusion exponent /r = 2. 
That estimate is obtained in a narrow region of the scal¬ 
ing variable, thus it may not be viewed as a true scaling 
exponent. However, taking it as an effective estimate 
of p, and taking a ~ 0.43 (from the scaling variable in 
Fig. 9), we obtain (z^) ^ [Eqs. ([T]) and ((2)]. This 
corresponds to an apparent subdiffusion. 

These results and those for walkers of type (I) in 
BDNNN deposits show important similarities. Both sys¬ 
tems have stretched exponential decay of the concentra¬ 
tion profile, a feature that is enhanced in the tube. They 
also show slowly increasing (subdiffusive) mean-square 
displacement, this effect being enhanced in the porous 
deposits. The slopes of scaling plots differ due to the 
different microscopic structures. 

Surprisingly, the effective exponents obtained from 
Fig. 9 (/r ~ 0.7, a ~ 0.3) are much closer to the ones 
for Pt atoms entering porous carbon samples (/r Ri 0.55 
and a ~ 0.2) Q, with the particularly interesting feature 
of highly stretched concentration profile. This suggests 
that the anomalous scaling observed in Pt atom diffusion 
may be a consequence of a very long crossover due to the 
confined geometry of those samples. 
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C. Tubes with exponentially decreasing cross 
section 

Here we consider the case of exponential decay ot A {z) 
because this facilitates the analytical study of crossover 
effects: 


A{z) = AQeyi^{-z/zo). (11) 

This gives N {z, t) = Ni exp [— 2 ^/ {2Dt) — z/zq] , where 
Ni depends on Aq, zq, and t should not be confused 
with Nq in Eq. 0: Nq = Ni exp {—z/zq)]. The crossover 
position is 


Zc ~ 2Dt/zo 


( 12 ) 


and the concentration relative to the origin at the 
crossover region is 


Pc = N {zc, t) /N (0, t) ~ exp {-ADt/zQ^). 


(13) 


For z ^ Zc-, the asymptotically normal diffusion is ob¬ 
served, with density p pc- This corresponds to very 
small concentrations at long times, thus this regime is 
difficult to be observed in simulation or in possible ex¬ 
perimental realizations. For z Zc, N {z,t) follows the 
simple exponential decay of A{z), which is the templat- 
ing effect of the tube shape. The diffusion seems to be 
anomalous, with p k, 1, and the concentration is much 
larger than pc- As time increases, this regime extends to 
larger distances and lower densities. 

The mean-square displacement (z^ (t)) ^ 

z^p{z,t)dz is estimated by finding the 
maximal value of 


z‘^p[z,t), which 


zo » 
(z2) . 


occurs for 

1 / Dt/zQ^ -I- 8 — '/dI/zq \/~Di/2. For 

y/m ~ Zdiff [Eq. (ITOl) ]. normal scaling with 
' t is obtained. This is a case in which the 
tube does not template the concentration profile. For 
Zo <C a/IV ^ Zdiff, (z^) ~ zo^, with corrections in 1/t, 
since the concentration profile follows the tube shape 
\p ^ A{z) in the region contributing to (z^)]. This is 
a case of anomalously slow diffusion (a = 0). These 
limiting behaviors indicate that subdiffusion is expected 
in the crossover region. 

The crossover is illustrated in Fig. 10a for a tube 
with Zo = 20 at three different times. The distributions 
were obtained in simulations of 10^ different walks with 
D = 1/3. For t = 10^, a crossover is observed from an 
initial regime with slope ~ 0.6 to a final regime with 
slope near 2. Again, we stress that these values are ef¬ 
fective exponents representative only of a narrow scaling 
region. Eq. (fT^ predicts the crossover length Zc « 33 
for Zo = 20 and t = 10^, in good agreement with the 
results in Fig. 10a. For t = 10®, Fig. 10a does not show 
any signature of a crossover up to z ~ 180, which is the 
maximal depth in which accurate concentrations could 
be estimated. Moreover, the data in Fig. 10a do not fit a 
single scaling curve with any variable of the form z/t“/^. 




ln(z/t°'‘") 


FIG. 10: (Color online) (a) Normalized walker concentra¬ 
tion as a function of the depth in a tube with exponentially 
decreasing A(z) with zq = 20 at t = 1000 (red squares), 
t = 10000 (green crosses), and t = 100000 (blue triangles). 
The dashed line has slope 0.6 and the solid line has slope 
2. (b) Normalized walker concentration as a function of the 
scaled depth in a tube with exponentially decreasing A{z) 
with Zo = 50 at t = 1000 (red squares), t = 2000 (green 
crosses), and t = 5000 (blue triangles). The dashed line has 
slope 0.6 and the solid line has slope 1.7. 


An illustration of the crossover scaling is possible in 
tubes with zq = 50 and a smaller time range. Fig. 
10b shows In [—In (A^/A/))] as a function of lnz/t°-^® at 
t = 1000, t = 5000, and t = 10000. The traditional form 
of the scaling variable gives a reasonable data collapse for 
large z in this system and indicates an effective anoma¬ 
lous diffusion (note that a complete scaling plot in this 
problem would have to involve the crossover length z^). 
The initial slope is near 0.6 and crosses over to a slope 
near 1.7 for the largest values of z. These values are close 
to the predictions of effective exponents 1 and 2 of the 
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crossover scaling; deviations are justified by the limited 
range of z in which the concentration is not negligible. 
The crossover positions predicted by Eq. (j9|) range from 
2 ;c Ri 13 {t = 1000) to Zc ~ 130 {t = 10000), which are 
close to the points with largest curvatures in the plots of 
Fig. 10b. 

V. CONCLUSION 

In the first part of this work, we studied the statistics of 
RW starting at the outer surface of porous deposits and 
restricted to move along their internal walls. BD and 
BDNNN models were used to produce the porous solids. 
These conditions are chosen to parallel the penetration 
of metal atoms in porous materials during the deposition 
of a film at their external surfaces. The concentration 
decay with the depth z is slower than Gaussian at short 
times. At long times, a change to normal diffusion is ob¬ 
served in most cases. However, the stretched exponential 
concentration profile remains in the most porous solids, 
produced by BDNNN, if the walker steps are restricted 
to nearest neighbors. This feature is correlated to a de¬ 
cay of the area explored by the RW approximately with 
the inverse z. However, RW of tracer particles left at 
various points of the solid do not show a significant de¬ 
pendence of the diffusion coefficient with z. This rules 
out the description of the anomalous diffusion by equa¬ 
tions with position-dependent coefficients . This 

shows that local features in space or time are not of main 
importance for the anomaly, but the global history in the 
penetration of the diffusing particles is the main ingredi¬ 
ent. Starting from surface deposition sites, not all inner 
sites at a given height are accessible in an equivalent way, 
some sites having a preferential path connecting them to 
the outer surface. 

In the second part of the work, we propose a model of 
RW confined to a tube of decreasing cross section A{z). 
We first consider an area which decreases slower than 


a Gaussian, similar to the decay obtained in one of the 
simulated experiment of the first part and account for a 
subdiffusion regime. Then considering the specific case 
of A(z) with a simple exponential decay, we predict an¬ 
alytically the crossover from a subdiffusion to a normal 
diffusion regime. The crossover position and crossover 
density show good agreement with numerical data. In 
both cases, the anomalous scaling is understood as a 
templating effect of A(z) on the concentration profile, 
which is measured in solid slices of constant area, but 
with walkers exploring only a fraction of this area. 

The above results suggest an alternative explanation 
to the stretched exponential concentration profiles and 
anomalous diffusion shown in recent works on the pen¬ 
etration of Pt atoms in disordered porous deposits [ 3 ], 
despite the difference in the exponent values. An interest¬ 
ing feature of our models is the independence of diffusion 
coefficients on the position, which is related to the lo¬ 
cal homogeneiW of the media, while the model equations 
used in Ref. [ 3 , 0 assume a dependence of those coef¬ 
ficients in time and position. On the other hand, those 
experiments also include a flux of atoms to the outer sur¬ 
face of the porous media, which is a feature not included 
in the present models and certainly important for their 
quantitative description. 

Independently of this possible application, the present 
work highlights the difficulty in interpreting experimental 
or computational data on anomalous diffusion, particu¬ 
larly if only the concentration profiles are measured or 
large crossover times are present. This is in agreement 
with recent experimental and theoretical works on the 
subject S ESI, HI. 
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